Opto- and electro-mechanical entanglement improved by modulation 
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One of the main milestones in the study of opto- and electro-mechanical systems is to certify entanglement 
between a mechanical resonator and an optical or microwave mode of a cavity field. In this work, we show 
how a suitable time-periodic modulation can help to achieve large degrees of entanglement, building upon the 
framework introduced in [Phys. Rev. Lett. 103, 213603 (2009)]. It is demonstrated that with suitable driving, the 
maximum degree of entanglement can be significantly enhanced, in a way exhibiting a non-trivial dependence 
on the specifics of the modulation. Such time-dependent driving might help experimentally achieving entangled 
mechanical systems also in situations when quantum correlations are otherwise suppressed by thermal noise. 



INTRODUCTION 

Opto-mechanical flJQ and electro-mechanical systems ||8]- 
[T3l are promising candidates for realizing architectures ex- 
hibiting quantum behavior in macroscopic structures. Once 
the quantum regime is reached, exciting applications in quan- 
tum technologies such as realizing precise force sensors are 
conceivable | IS TSl. One of the requirements to render such 
an approach feasible, needless to say, is to be able to certify 
that a mechanical degree of freedom is deeply in the quantum 
regime ll6H20il . The detection of entanglement arguably con- 
stitutes the ultimate benchmark in this respect. While effective 
ground state cooling has indeed been experimentally closely 
approached ||6l[T0l and achieved QSIII], the detection of 
entanglement is still awaiting. 

In this work, we emphasize that a mere suitable time- 
modulation of the driving field may significantly help to 
achieve entanglement between a mechanical mode and a ra- 
diation mode of the system. We extend the idea of Ref. 
jJTl, putting emphasis on the improvement of entanglement 
by means of suitable modulations |21-23|. The method used 
here is not a direct modulation of the frequencies of the two 
modes (parametric amplification), but the system is instead 
externally driven with a modulated field. This time depen- 
dence of the driving indirectly affects the effective radiation 
pressure coupling between the two modes and generates non- 
trivial entanglement resonances. In this way, with the appro- 
priate choice of the modulation pattern, large degrees of two- 
mode squeezing can be reached. 



MODULATED OPTO- AND ELECTRO-MECHANICAL 
SYSTEMS 

We consider the simplest scenario of a mechanical res- 
onator of frequency a;^ coupled to a single mode of the elec- 
tromagnetic field of frequency uja- This radiation field could 
be an optical mode of a Fabry-Perot cavity [T8l[T9ll24l 
or a microwave mode of a superconductive circuit ll8 molfT4l . 
It can be shown that the Hamiltonians associated to this two 
experimental settings are formally equivalent lfT4l [T9l and 
therefore the theory that we are going to introduce is general 



enough to describe both types of systems. 

We assume that the radiation mode is driven by a coherent 
field with a time dependent amplitude E{t) and frequency a;;. 
The particular choice of the time dependence is left unspeci- 
fied but we impose the structure of a periodic modulation such 
that E{t + t) = E{t) for some r > of the order of uj:^^. 
In this sense, the driving regime that we are going to study is 
intermediate between the two opposite extremes of constant 
amplitude and short pulses. The Hamiltonian of the system is 

H = hhJao) a + -hojmip^ + Q^) — figo) aq 
+ih[E{t)e-"^'*'a^ - E*{t)e"^'^a], (1) 

where the mechanical mode is described in terms of dimen- 
sionless position and momentum operators satisfying [q^p] = 
i, while the radiation mode is captured by creation and an- 
nihilation operators obeying the bosonic commutation rule 
[a, a^] — 1. The two modes interact via a radiation pressure 
potential with a strength given by the coupling parameter g. 

In addition to this coherent dynamics, the mechanical mode 
will be unavoidably damped at a rate 7™, while the opti- 
cal/microwave mode will decay at a rate k. These dissipative 
processes and the associated fluctuations can be taken into ac- 
count in the Heisenberg picture by the following set of quan- 
tum Langevin equations lfT4l[T7l - [l9l . 

q = ujjnP, (2) 

p = -UJrnq - IniP + 90^0- + ^, 

a = - (k + iA)a + igaq + E{t) + \/2Ka". 

In this set of equations a convenient rotating frame has be 
chosen a i~> ae^*"'*, such that the detuning parameter is 
A = LUa — ijOi- The operators ^ and a'" represent the me- 
chanical and optical bath operators respectively, and their cor- 
relation functions are well approximated by delta functions 

{mi{t') + e(t')e(t)>/2 = 7m(2n™ + l)5{t - t'), (3) 
{c^^{ty^\t'))^{n,^l)5{t~t'), 

where rix — {exp{hijjx / {ksT)) — 1)^^, is the bosonic mean 
occupation number at temperature T. 
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CLASSICAL PERIODIC ORBITS: FIRST MOMENTS 

We are interested in the coherent strong driving regime 
when (a) ^ 1. In this limit, the semiclassical approximations 
(a^a) ~ and (aq) ~ («)('?) are good approximations. 

Within this approximation, one can average both sides of Eq. 
Q and get a differential equation for the first moments of the 
canonical coordinates 



{p) = -^m{q) - im{p) + .g|(a)P, 

(a) = -(K + iA)(a) +ig(a){q) + E{t). 



(4) 



Far away from the well known opto- and electro- 
mechanical instabilities, asymptotic r-periodic solutions can 
be used as ansatz for Eqs. Q (see the Appendix for a more 
detailed analysis). These solutions represent periodic orbits in 
phase space and are usually called limit cycles. These cycles 
are induced by the modulation and should not be confused 
with the limit cycles emerging in the strong driving regime 
due to the non-linearity of the system. Because of the asymp- 
totic periodicity of the solutions, one can define the funda- 
mental modulation frequency as = 2tt/t, such that each 
periodic solution can be expanded in the following Fourier se- 



(0(t))- ^ 0„e""*, = q,p,a. 



(5) 



The Fourier coefficients {0„} appearing in Eq. (jsj) can be ana- 
lytically estimated as shown in Appendix and they completely 
characterize the classical asymptotic dynamics of the system. 

Finally we notice that the classical evolution of the dynam- 
ical variables will shift the detuning to the effective value of 
A(f) = A—g{q{t)). For the same reason, it is also convenient 
to introduce an effective coupling constant defined as 



Since we are in the limit in which classical orbits emerge 
((a) ^ 1), it is a reasonable approximation to express the 
previous set of Langevin equations Q in terms of the new 
fluctuation operators (j7]i and neglect all quadratic powers of 
them. The resulting linearized system can be written as a ma- 
trix equation [21 1, 



where. 



A{t) 



u = A(t)u + n{t), 





-l^m -7m '^git) '^Kt) 

-^g{t) -K A{t) 

^g{t) -A{t) -K 



(9) 



(10) 



is a real time-dependent matrix. 

If the system is stable, and as long as the linearization is 
valid, the quantum state of the system will converge to a Gaus- 
sian state with time dependent first and second moments. The 
first moments of the state correspond to the classical limit cy- 
cles introduced in the previous section. The second moments 
can be expressed in terms of the covariance matrix V{t) with 
entries 

Vk,i{t) = {uk{t)u\{t) + u\{t)uk(t)) /2. (11) 
One can also define a diffusion matrix D as 

6{t - t')Dkj = {nk{t)n\{t') + n|(i')nfc(<))/2, (12) 

which, from the properties of the bath operators (|3]l, is diago- 
nal and equal to 

D = diag[0, 7(2n„ + 1), K{2na + 1), «;(2n„ + 1)]. (13) 



m = wHt))/^. 



(6) 



From Eqs. (|9]l and ( 12 1, one can easily derive a linear differ- 
ential equation for the correlation matrix. 



QUANTUM CORRELATIONS: SECOND MOMENTS 

The classical limit cycles are given by the asymptotic solu- 
tions of Eqs. (|4]). In order to capture the quantum fluctuations 
around the classical orbits, we introduce a column vector of 
new quadrature operators u — [5q, Sp, 6x, SyY" defined as: 



Sq = q- {q{t)), 
Sp = p- {p{t)), 

6x^ [{a - (ait))) + {a - {a{t)))^ /V2, 
5y^-z[ia-{am-ia-{a{t))y] /^/2. 



(7) 



This set of canonical coordinates can be viewed as describing 
a time-dependent reference frame co-moving with the classi- 
cal orbits. The corresponding vector of noise operators will 
be 



[0,e,(a" 



2, -i{a" 



(8) 



-Vit) = Ait)Vit) + Vit)A'^{t) + D. 



(14) 



Since the first and the second moments are specified, Eqs. Q 
and ( [T4| provide a complete description of the asymptotic dy- 
namics of the system. Apart from the linearization around 
classical cycles, no further approximation has been done: Nei- 
ther a weak coupling, adiabatic or rotating-wave approxima- 
tion. Numerical solutions of both equations (j4]) and ( 14 1 can 
be straightforwardly found. These solutions will be used to 
calculate the exact amount opto- and electro-mechanical en- 
tanglement present in the system. 

The asymptotic periodicity of the classical solutions (Eq. 
^) implies that, in the long time limit, A{t + T) = A{t). This 
means that Eq. ([T4]l is a linear differential equation with peri- 
odic coefficients and then all the machinery of Floquet theory 
is in principle applicable. Here, however, since we are only 
interested on asymptotic solutions, we are not going to study 
all the Floquet exponents of the system. The only property 
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that we need is that, in the long time limit, stable solutions 
will acquire the same periodicity of the coefficients: 



V{t + T) = V{t). 



(15) 



This is a simple corollary of Floquet's theorem. In the sub- 
sequent sections we will apply the previous theory to some 
particular experimental setting and show how a simple modu- 
lation of the driving field can significantly improve the amount 
of opto- and electro-mechanical entanglement. 



ENTANGLEMENT RESONANCES 

In this section we are going to study what kind of amplitude 
modulation is optimal for generating entanglement between 
the radiation and mechanical modes. As a measure of entan- 
glement we use the logarithmic negativity En which, since 
the state is Gaussian, can be easily computed directly from 
the correlation matrix V{t) [26-28 |. We have also seen that 
the correlation matrix is, in the long time limit, r-periodic. 
This suggests that it is sufficient to study the variation of en- 
tanglement in a finite interval of time [t, t + t] for large times 
t. One can then define the maximum amount of achievable 
entanglement as 



E 



N 



lim max Ef^[h). 

t-s-oo h£[t,t+T] 



(16) 



This will be the quantity that we are going to optimize. 

We first study a very simple set of parameters (see caption 
of Figure [TJ in order to understand what the optimal choice is 
for the modulation frequency. For this purpose we impose the 
effective coupling to have this simple structure 



9{i) = .90 + gn 



-iflt 



(17) 



where go is associated to the main driving field with detun- 
ing A, while gQ is the amplitude of a further sideband shifted 
by a frequency from the main carrier. Without loss of gen- 
erality we will assume go and 50 to be positive reals. This 
kind of driving is a natural one and has been chosen for rea- 
sons that will become clear later From now on we set the 
detuning of the carrier frequency to be equal to the mechan- 
ical frequency A — uj,n. This choice of the detuning cor- 
responds to the well known sideband cooling setting ifTTl l24l 
and it has been shown to be also optimal for maximizing opto- 
mechanical entanglement with a non-modulated driving lfT9l . 
Fig.[l]shows the maximum entanglement Ej^ between the me- 
chanical and the radiation modes as a function of the modu- 
lation frequency O and for different values of the driving am- 
plitude go. This maximum degree of entanglement has been 
calculated for t > 200/k when the system has well reached 
its periodic steady state. 

We observe that in Fig.[T|there are two main resonant peaks 
at the modulation frequencies 
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FIG. 1. Maximum entanglement En as a function of the modulation 
frequency Q and for different values of the driving strength go . The 
chosen parameters in units of cjm are: k — 0.2, 7™ = 10"", A = 1, 
rim = na = 0, gn = 0.1, go = 0.2 (circles), 0.4 (squares), 0.6 
(diamonds), 0.8 (triangles). 



We will now provide some intuition why one should expect 
the main resonances at the locations where they are observed. 
First assume that go — 0. Then, for A = Wm, the linearized 
Hamiltonian in the interaction picture is 

{e^^'-Sb'' +Sbe'"^"^) /2, (19) 

where the bosonic operators are defined as 5a — {5x + 
iSy)/V^, Sb = {Sq + i6p)/y/2. From Eq. (19i, it is clear 
that for n = 2a;„i, neglecting all rotating terms, we get the 
well known two-mode squeezing generator 



int 



-hgn (JaUfe^ + SaSb) /2. 



(20) 



So, in the case of go = 0, a modulation of il = 2uj„i would be 
the most reasonable choice in order to generate entanglement. 
However, this regime is well known to be highly unstable and, 
in practice, it cannot be used for preparing entangled steady 
states (14]. 

This is why we need to consider a modulated coupling of 
the form given in Eq. ( 17 1 - or a similar type of modulation 



sharing these features. We now allow for go being different 
from zero, giving rise to a situation which can be assessed in 
a very similar way as above (only that the rotation terms will 
take a more involved form). The main amplitude go then takes 
the role of cooling and stabilizing the system while the mod- 
ulation amplitude gn is used to generate entanglement. At the 
same time however, as shown in Refs. ||5]|25l, for go > k/-\/2 
the system hybridizes in two normal modes of frequencies 



± .90/2. 



(21) 



As a consequence, this will affect the modulation frequency 
n that one has to choose in order to achieve the two-mode 



n ~ 2u!jn ± 50- 



(18) 



squeezing interaction given in Eq. (20i. This is the reason 
for the presence of two resonant peaks in Fig. [T] and for the 
resonance condition given in Eq. ( 18 1. 
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Note also that the choices of modulations that give rise 
to the optimal local single-mode squeezing 121] of the me- 
chanical mode and the degree of entanglement are not iden- 
tical. This is rooted in the "monogamous nature" of squeez- 
ing: For a fixed spectrum of the covariance matrix, one can 
either have large local or two-mode squeezing. This effect is 
observed when considering the modulation frequencies that 
achieve maximum single- and two-mode squeezing. 

We finally observe that the height of the two peaks, due 
to the cavity filtering, is not equal: the first resonance at 

= 2Ldm — go is better for the amount of steady state en- 
tanglement. One could also ask what the behavior of entan- 
glement is when we change the amplitude of the modulation. 
Fig. [2] shows the amount of entanglement En as a function 
of cjQ and for different choices of go. We observe that entan- 
glement is monotonically increasing in g^ up to a threshold 
where the system becomes unstable. 
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FIG. 2. Maximum entanglement En as a function of the modulation 
amplitude gn and for different values of the driving strength ^0. The 
chosen parameters in units of Wm are: k = 0.2, 7™ = 10~®, A = 1, 
rim — ria = 0, — 2ujm — go, go = 0.2 (circles), 0.4 (squares), 
0.6 (diamonds), 0.8 (triangles). 



OPTO- AND ELECTRO-MECHANICAL ENTANGLEMENT 
IN REALISTIC SETTINGS 

We have seen that an effective coupling of the form g{t) = 



go + gn e 



-i(2a;„j-go)t 



is optimal for the generation of entan- 



glement within the considered class of drivings. However, the 
parameter g{t) depends on the average amplitude (a(t)) and 
assuming such a simple structure may seem somewhat artifi- 
cial. In this section, we show how the desired time-dependent 
coupling can indirectly result from the classical limit cycles of 
the system (see insets of Figs. |3] and |4|i and we also take into 
account the effect of a temperature of the order of T ~ 100 
mK. The natural "educated guess" for the structure of the driv- 
ing field will be 



E{t) ^Eo + EnEe 



-i{2uim-go)t 



(22) 



For the choice of the other parameters, we focus on two set 
of parameters corresponding to two completely different sys- 



tems: an optical cavity with a moving mirror and a supercon- 
ducting wave guide coupled to a mechanical resonator The 
parameters are chosen according to realistic experimental set- 
tings, see, e.g., Ref. |5 1 (opto-mechanical system) and Ref. ||9| 
(electro-mechanical system). Fig. [3] and Fig. |4] show that, in 
both experimental scenarios, entanglement can significantly 
be increased by an appropriate modulation of the driving field. 
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FIG. 3. (Optical cavity). The degree of entanglement, measured in 
terms of the logarithmic negativity, as a function of time. The full line 
refers to a modulated driving (fl = lAum) while the dotted line cor- 
responds to a non-modulated driving (SI = 0). The chosen parame- 
ters in units of cjrn are: k — 0.2, jm = lO"'^, A = 1, rim — 2x 10'', 
ria = 0, go = 4 X 10~^ Eo = 7 x lO"*, En = 2.5 x 10*. The inset 
shows the trajectory of the effective coupling g{t) = \/2g{a{t)) in 
the complex plane due to the time evolution of the optical amplitude. 
The phase space orbit (black line) is numerically simulated from Eq. 

while the limit cycle (green line) is an analytical approximation 
(see Appendix for more details). 
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FIG. 4. (Microwave cavity). Entanglement log-negativity as a func- 
tion of time. The full line refers to a modulated driving (fi = 
l.Sum) while the dotted line corresponds to a non-modulated driv- 
ing (Jl = 0). The chosen parameters in units of ujm are: k — 0.02, 
7m = 3 X 10~^ A = 1, rim = 200, ria = 0.03, 50 = 2 x 10"^ 
So = 9 X 10^, En = 1.3 x 10^. The inset depicts the trajectory of 
the effective coupling g{t) — \f2g{a{t)) in the complex plane due 
to the time evolution of the microwave amplitude. The phase space 
orbit (black line) is numerically simulated from Eq. (j4]l, while the 
limit cycle (green line) is an analytical approximation (see Appendix 
for more details). 
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SUMMARY 

In this work, we have shown how time-modulation can 
significantly enhance the maximum degree of entanglement. 
Triggered by the time-modulated driving, the mode of the 
electromechanical field as well as the mechanical mode start 
"rotating around each other" in a complex fashion, giving rise 
to increased degrees of entanglement. The dependence on the 
frequencies of the additional modulation is intricate, with res- 
onances highly improving the amount of entanglement that 
can be reached. The ideas presented here could be particu- 
larly beneficial to prepare systems in entangled states in the 
first place, in scenarios where the parameters are such that the 
states prepared are close to the boundary to entangled states, 
but where this boundary is otherwise not yet quite reachable 
with present technology. At the same time, such ideas are 
expected to be useful in metrological applications whenever 
high degrees of entanglement are needed. 



of recursive relations: 



<Zr!,0 — Pn,0 — 0, a„_o — 



K + i(Ao + nil) ' 



(26) 



corresponding to the 0-order perturbation with respect to Go, 
and 



Pn,j — Qn , J : 



(27) 



j — 1 OO ^ 

k—Q m — — OO '"' ' 

j-1 OO 

. \ ^ \ ^ ^m,kQ.n—in.j — k~l 



(29) 



giving all the j-order coefficients in a recursive way. For all 
the examples analyzed in this paper we truncated the analyti- 
cal solutions up to j < 3 and \n\ < 2. This level of approx- 
imation is already high enough to well reproduce the exact 
numerical solutions. 



APPENDIX 

In this appendix we derive analytical formulas for the 
asymptotic solutions of the classical system of dynamical 
equations Q. A crucial assumption for the following pro- 
cedure is that it is possible to expand the solutions in powers 
of the the coupling constant 50 
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(23) 



3=0 



where O — a,p,q. This is justified only if the system is 
far away from multi-stabilities and the radiation pressure cou- 
pling can be treated in a perturbative way. A very important 
feature of the set of equations Q is that they contain only two 
non linear terms and those terms are proportional to the cou- 
pling parameter go- This implies that, if we use the ansatz 



(23 I, each function Oj will be a solution of linear differen- 
tial equation with time dependent parameters depending on 
the previous solution Oj_i(<). Since E{t) = E{t + t), from 
a recursive application of Floquet's theorem, follows that sta- 
ble solutions will converge to periodic limit cycles having the 
same periodicity of the driving: {0{t)) — {0{t + r)). One 
can exploit this property and perform a double expansion in 
powers of go and in terms of Fourier components 



00 00 



(24) 



j—O n— — 00 



where n are integers and ft — 2tt/t. A similar Fourier series 
can be written for the periodic driving field, 



(25) 



n— — 00 



The coefficients Onj can be found by direct substitution in 
Eq. Q. They are completely determined by the following set 
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